Interferometric detection of a single vortex in a dilute Bose-Einstein condensate 
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Using two radio frequency pulses separated in time we perform an amplitude division interference 
experiment on a rubidium Bose-Einstein condensate. The presence of a quantized vortex, which is 
nucleated by stirring the condensate with a laser beam, is revealed by a dislocation in the fringe 
pattern. 



The experimental realization of Bose Einstein conden- 
sation of dilute atomic gases has opened for study a whole 
new class of macroscopic quantum systems . In partic- 
ular, the question of the superfluidity of a dilute conden- 
sate can now be addressed experimentally by analyzing 
the response of such a system to a moving perturbation. 
This response is severely restricted by the strong con- 
straints that exist on the velocity field of the condensate. 
For a macroscopic wave function tp{r) = \/n{r) e'^^^^^ 
where n(r) is the density of the condensate and S{r) the 
local phase, the velocity field is {h/M)V S ^ where M is 
the mass of a particle of the fluid. This relation requires 
that the circulation of the velocity field around any closed 
contour be quantized and equal to nh/M^ where n is an 
integer . A non-zero value of n is associated with the 
presence of one or several quantum vortices. Such a vor- 
tex is characterized by a line along which the condensate 
density vanishes and around which the phase of the wave 
function varies by 2ti7i (for recent reviews, see [||J|]). 

A reliable way to nucleate a quantum vortex in a di- 
lute Bose-Einstein condensate is to place the system in 
a rotating potential as was originally done with super- 
fluid helium in the famous "rotating bucket experiment" 
Starting with a condensate confined in a cylindri- 
cally symmetric trap, one can stir the atomic cloud us- 
ing an additional rotating anisotropic potential created 
by the dipole potential of a non-resonant stirring laser 
beam [||J|] . There exists a range of values of the stirring 
frequency for which one can nucleate in this way a single 
vortex which will center itself on the symmetry axis of 
the trap. The presence of the vortex is then revealed by 
the existence of a density dip at the center of the spatial 
distribution of the condensate. Further experiments have 
also revealed that the angular momentum per particle is 
equal to % for this centered, one- vortex state |]lQ| , pj] . 

We present in this article a direct observation of the 
phase of such a rotating condensate. We follow the 
scheme proposed in which relies on an interfer- 

ometric measurement. Here, we divide the wave func- 
tion of the condensate into two components separated in 
position and momentum. In the location where these 
components overlap, a matter- wave interference fringe 
pattern is produced which reveals the phase difference 



between the two components. This homodyne detection 
of the phase pattern of the vortex state of the conden- 
sate is analogous to the experiments aiming at studying 
the transverse phase structure of optical Laguerre-Gauss 
modes jl^. Moreover, it is complementary to the het- 
erodyne detection of a vortex state performed in Boulder 
in 1999 In that experiment, a vortex was created in 
a two component condensate by a direct printing of the 
vortex state phase e*^ on one of the two components {0 
is the azimuthal angle around the symmetry axis of the 
condensate). The second component, corresponding to a 
different internal state of the atoms, remained at rest at 
the center of the confining potential and was used as a 
"reference" matter wave to detect the phase of the vortex 
state contained in the rotating component. In this case, 
the two interferring components correspond to two am- 
plitudes with a different spatial phase while in our case 
the two components have the same spatial phase but are 
shifted in position and momentum. 

We prepare our Bose-Einstein condensate by evapora- 
tive cooling of 10^ rubidium (^^Rb) atoms confined in 
a loffe-Prit chard trap. The trap is axisymmetric and 
its eigenfrequencies are respectively uo±_ = 27r x 192 Hz 
and uoz = 27r x 11.7 Hz for the stretched state m = +2 
of the 5S1/27F = 2 ground state. The initial tempera- 
ture of the atom cloud, which is precooled using optical 
molasses, is 100 /iK. At the end of the evaporation, a 
condensate containing 3 10^ atoms is formed, and the 
temperature of the cloud is below 100 nK. We then nu- 
cleate the vortices by the method described in ||] : we su- 
perimpose to the isotropic magnetic trap a far detuned 
laser beam propagating along the weak horizontal axis 
{z) of the trap and positioned by two crossed acousto- 
optic modulators. This light beam creates in the xy 
plane an anisotropic optical dipole potential of the form 
5U{X,Y) =eMuo\{X'^ - Y'^)/2. We rotate at constant 
frequency ^ 130 Hz the eigenaxes XY of this po- 

tential, which provides the desired stirring, with e = 0.06 
in the present experiment. We stir the atomic cloud for 
a duration of 300 ms. We verified in a preliminary step 
that this duration is sufficiently long to nucleate a well 
centered single vortex. 

To realize the amplitude splitting interferometer, we 
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use a method similar to the one described in [T^]. A 
succession of two radio frequency (rf) pulses, which reso- 
nantly drive transitions between different internal states 
of each rubidium atom, is used to coherently split and 
recombine the condensate cloud. The atomic rubidium 
condensate is initially polarized the F = m = 2 Zeeman 
substate. The first rf pulse couples this magnetic spin 
state to the four other Zeeman substates of the F = 2 
manifold. Among these four, only the m = +1 is also 
magnetically trapped. However, due to gravity, the cen- 
ter of its confining potential is vertically displaced with 
respect to the center of the m = +2 trap by an amount 
Axq = g/uo\^ Q fim. After this first pulse, the atoms are 
allowed to evolve in the harmonic trapping potential for 
a variable time ri. During this time, those atoms driven 
into the other Zeeman substates fall or are ejected out 
of the trap. What remains are the two trapped clouds 
now with slightly different momenta and positions as de- 
termined by their respective phase space trajectories pic- 
tured in Fig. I^a. When the second rf pulse is applied, 
each trapped state generates a copy in the other trapped 
Zeeman sublevel making a total of four wave packets 
present in the trap (two in m = +1 and two in m = +2). 
We then let the atoms evolve in the harmonic poten- 
tial for a time r2 (Fig. [l]b) after which we extinguish the 
trapping potential and let the condensate clouds freely 
expand and fall for 25 ms. Finally, we record the atomic 
interference pattern between the overlaping wave packets 
of a given internal Zeeman spin state using the resonant 
absorption imaging technique originally used to detect 
the density dip associated with the vortex line. 




FIG. 1. Phase-space trajectories of the m = +2/m = +1 
Zeeman substate wave packets in the harmonic trapping po- 
tential during the experimental sequence, (a): location of the 
two wave packets at a time n after the first rf pulse, (b): lo- 
cation of the four confined wave packets at a time T2 after the 
second rf pulse. The hatched disks correspond to the wave 
packets generated by the second rf pulse. 

The possibility of changing independently n and r2 
allows for an independent control of the initial relative 
velocity and separation of two interfering wave packets. 
Note that the amplitude of these rf pulses is large enough 
so as to resonantly excite all atoms in the condensate, in 



spite of the inhomogeneity of the magnetic field inside 
the trap. This requires that the Rabi frequency Q^f asso- 
ciated with these rf excitations be at least of the order of 
/jj/h, where fi is the chemical potential. In our case, fi/h 
is typically 27r x 10 kHz while ftrf = 27r x 7 kHz. The 
duration of each rf pulse is 13 jas which ensures that the 
populations of the two states m = +2 and m = +1 after 
each pulse are comparable. 

Before showing the interference patterns detected ex- 
perimentally, we present a simple theoretical model to 
discuss the type of signal which can be expected. Con- 
sider first a trapped condensate localized around r = 
with no velocity field ('^(r) = yn(r)). We assume 
that this condensate is split into two parts separated 
by a, so that the unnormalized state of the system is 
?/^(r— a/2)+?/^(r+a/2). The potential confining the atoms 
is then extinguished and the two clouds overlap and inter- 
fere. The detection of the interference pattern is done at 
a time t large enough so that the final size of the conden- 
sate is much larger than the initial one. Consequently a 
model relying on a non-interacting Bose-Einstein conden- 
sate is sufficient to interprete the fringe pattern [^9|-^. 
The phase pattern of the wavefunction issued from the 
condensate initially localized in tq = ±a/2 is the same 
as that of the function 

^(r)=exp^^^j (1) 

so that the interference pattern consists in parallel 
fringes, with a fringe spacing p2| : 

The phase pattern of (|l|) remains valid even if the cloud 
initially localized in tq has a non zero mean velocity. 
Therefore the expression for the fringe spacing (1) is still 
correct when the relative velocity of the two interfering 
clouds is not zero, provided a denotes the initial spatial 
separation between the centers of the condensates before 
the atoms were released [^. Note that the size of the 
expanding condensates also increases linearly with t, so 
that the number of visible fringes is independent of the 
expansion time. In practice, the number of visible fringes 
is limited by the optical resolution Sx of the imaging sys- 
tem {Sx = 7 /im in our experiment), since one must have 
6x <^ Xs. Choosing the maximal value of the expansion 
time compatible with our experimental setup (25 ms) 
gives from the product x^a ~ (11 /im)^. Then for 
a good fringe visibility, we choose a small value of a, 
typically a = 2.5 /im, which yields Xs ~ 45 /im. The 
diameter of the condensate after expansion is 100 /im, 
so that 3 bright fringes should be clearly visible in the 
interference pattern. The corresponding calculated pat- 
tern is showed in Fig. |^a assuming an initial Gaussian 
distribution with an r.m.s. width Ax = Ay = 0.5 /im, a 
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separation a = 2.9 /im and an expansion time of 25 ms, 
yielding Xg = 39 /im. 

Consider now a condensate centered at r = tq, and as- 
sume that a quantized vortex, with an angular momen- 
tum -\-h with respect to the z axis, sits at the center of 
the condensate. The initial wave function of the conden- 
sate is ?/;(r) = {x-XQ^i{y-yQ)) f{\p-po\,z) e'*^""". Here 
we use cylindrical coordinates (p, z) with respect to the 
z axis, and / is a function of \p\ and z with a constant 
phase. When the condensate is released from the trap 
and expands for a time t, the phase pattern of the cloud 
wavefunction is equal to that of the function defined 
as: 

$'(r) = {x- x,+i{y - yi)) exp ( '^^^^^/°^' ) (3) 

where we set ri = tq + vgt. Note that we allow for a 
possible non zero velocity vq = h'k/M of the center of 
mass of the condensate. This result can be shown first 
for a cloud initially at rest (vq = 0), and can then be 
generalized to a moving cloud using a Galilean transform. 




FIG. 2. Expected fringe pattern of a Bose-Einstein con- 
densate initially splitted into two parts and undergoing a free 
expansion phase. Figure (a) is without a vortex and (b,c) are 
with a vortex. For (6), close to our experimental conditions, 
the fringe spacing Xs is equal to 39 /xm, and it is equal to the 
separation of the vortex cores after expansion |ri — r^i|. (c): 
same as (b), with a fringe spacing Xs = 13 /im= |ri — ri\/3 
(this fringe spacing is too small to be detected in our exper- 
imental setup). For figures (b) and (c), the relative phase of 
the two condensates is tt. 

The expected signal in an interference experiment can 
be obtained by superposing two wavefunctions varying 
as (||), obtained by expanding two condensates initially 
located at tq = a/2 and Tq = —a/2, with initial velocities 
Vq and Vq. The interference pattern consists in parallel 
lines separated by Xs = ht/Ma^ with "defects" associated 
with the presence of the vortex. These defects are clearly 
visible if the distance ri — r'l = a + (vq — v'o)^ between 
the two vortex cores after expansion is larger than Xg. 
Since a is very small compared to Xg in our setup, this 
visibility condition requires that: 



where Av = |vo — v'o|- On the other hand A'^ should 
be smaller than the velocity width of each condensate 



after expansion, so that the two clouds overlap after the 
time-of- flight. In practice, we choose the times ri and r2 
so that (Q) is close to being an equality and the overlap 
of the two wave packets after expansion is maximized. 
The centers of the two vortices after expansion are then 
separated by the fringe spacing Xs and one expects the 
characteristic fringe pattern shown in Fig. |^b, with two 
dislocations ( "H" shape) indicating the centers of the two 
vortices after time of flight. In Fig. |^, we show the 
fringe pattern which would be expected for a larger initial 
separation a, where more fringes are visible. 

We now turn to the discussion of the experimental re- 
sults. We have first measured the interference pattern in 
the absence of any vortices. As expected, this pattern 
consists in straight fringes, and a typical result is shown 
on Fig. ||. The two patterns corresponding to the m = +2 
and m = +1 components are recorded simultaneously. 
The spatial separation between these two components 
results from the difference between the initial velocities 
of the m = +1 and m = +2 clouds when released from 
the trap. The difference in the fringe spacings reflects 
the fact that the initial separation ai between the two 
wavepackets in state m = +1 was larger than the separa- 
tion a2 between the two wavepackets in m = +2 (see e.g. 
Fig. |l|). We have checked for each component that the 
variation of the fringe spacing with ai,2 is in good agree- 
ment with (H) (see also [Q). In the following we use 
the interference pattern obtained with the m = 2 com- 
ponent, which offers a better visibility due to its larger 
spatial extent after time of flight, which is itself due to a 
stronger confinement in the trap. 



m = +2 m = +l 



FIG. 3. Interference pattern of a non rotating condensate 
split into 4 parts (2 in each Zeeman state m = +2 and 
m = 1), which overlap during the free expansion phase. Here 
n — 0.3 ms and T2 — 0.4 ms. 

When we stir the atomic cloud at a frequency of 130 Hz 
(close to the critical frequency uj_\_/\^ l^), a well cen- 
tered, single vortex is created in the condensate. The 
result from an interference experiment performed in this 
case is shown on Fig. ^ for the m = +2 channel. This 
pattern, which is very different from the one obtained 
with a slightly slower rotation (125 Hz) at which no vor- 
tex nucleation occurs (Fig. ^), is the main result of this 
paper. It clearly shows the telltale "H" shape character- 
istic of the 27r phase winding due to the presence of a 
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vortex (see Fig. ^3). We repeated this experiment many 
times for the same initial condition, and have found that 
the interference patterns produced have a similar shape 
but with a relative phase between the two components 
which fluctuates from one shot to another so that the 
bright and dark regions in Fig. ^ can be exchanged (see 
e.g. Fig. ^). Finally, when several vortices are nucleated 
by stirring at a higher frequency, the fringe structure is 
more complex and the number of fringe dislocations is 
larger (Fig. ||d). 




FIG. 4. Interference pattern measured in the m=H-2 chan- 
nel with no (a), one (b-c) and several (d) vortices. For these 
pictures, n = 0.688 ms and T2 = 1.320 ms. The stirring fre- 
quency was set to Q = 27r X 125 Hz (a), Q = 27r x 130 Hz 
(b-c) and — 2iix 154 Hz (d). The patterns (b) and (c) were 
recorded with the same initial conditions, and the change in 
the interference pattern results from a change in the relative 
phase of the two parts of the condensate. 

In conclusion, we have observed the 27r phase shift as- 
sociated with the presence of a vortex in a single com- 
ponent Bose-Einstein condensate by using an interfero- 
metric detection method. This technique can be used as 
a very sensitive tool to diagnose the presence of vortices 
generated by the fast motion of a macroscopic object in 
a condensate. This was demonstrated very recently at 
MIT, in an experiment where the interference pattern 
between a condensate with vortices and a condensate at 
rest has been recorded [Q. 

The "beam splitter" used here is very simple, since it 
consists of a resonant radio frequency pulse. However, 
one disadvantage of this type of "beam splitter" is the 
unavoidable loss of atoms to the other untrapped Zee- 
man states resulting in at most only ~ 1/3 of the atoms 
participating in the interference signal in the m = +2 
channel. We plan to use in the near future a more ef- 



ficient scheme discussed in Q, which involves the use 
of two Bragg pulses to couple just two states of different 
momentum. 
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